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We derive radiation reaction directly from the S-matrix of QED. We identify the diagrams and 
processes contributing to recoil effects in the average momentum of a scattered electron and, in 
an explicit example, find exact agreement with classical electrodynamics in the limit fi — > 0. Our 
results unify a number of disparate, previously reported approaches. 



Understanding radiation reaction presents one of the 
oldest problems in electrodynamics, and has recently seen 
a renewal of interest due to the potential impact, and 
detection, of recoil effects in high-intensity laser-matter 
interactions [2 [2] . The years have seen many proposals 
for equations which describe a classical, radiating par- 
ticle and avoid the problems of the Lorentz-Abraham- 
Dirac ('LAD') equation [MS], but without consensus [7]. 
Nor is it completely understood how radiation reaction 
('RR') emerges from QED Q^. For example, the pho- 
ton spectrum emitted from an electron in a background 
field gives, in the classical limit, the radiation spectrum 
of a particle moving according to the Lorentz equation; 
no recoil effects are seen [H [9] . This has prompted pro- 
posals for which Feynman diagrams contain RR [TU , and 
the claim that quantum and classical electrodynamics are 
not compatible [llj . 

Given the great interest in this topic, we will here de- 
rive both quantum and classical RR directly and unam- 
biguously from QED. In doing so we will unify several 
existing approaches (such as use of the in-in formalism 
and coherent vs incoherent emissions), and resolve the 
above issue with the photon spectrum. Rather than pur- 
sue an equation, see e.g. [T7 for quantum stochastic cor- 
rections to LAD, wc will consider scattering in QED, ex- 
plain which diagrams contribute to RR, and illustrate 
these statements with an explicit example. All we will 
do is calculate (the expectation value of) the momen- 
tum of a quantum particle, and take h to zero to obtain 
the classical momentum. For our example, we consider 
scattering in a plane wave background modelling, say, an 
intense laser. With this choice, the Lorentz force can be 
treated without approximation, allowing a clean separa- 
tion of recoil (the impact of radiation on the electron's 
motion) and non-recoil effects. 

This paper is organised as follows. We begin by review- 
ing how and where RR arises classically. We then identify 
where RR arises in QED, stating which diagrams con- 
tribute. The electron recoil is then written down, and we 
take /i — 0, finding agreement with classical predictions. 
Connections to other approaches are then discussed be- 
fore we conclude. 

Classical:- The LAD equation for a radiating particle 
with orbit x''(t), proper time t, is {c = = 1, h ^ 1) 

2 

mXf, = eF^f{x)x'' + - — (x^i^, - i^x^)i'^ , (1) 



where F°^^ is an external field jj^. Retaining only the 
first term on the right hand side of ([T]) gives the Lorentz 
force equation, the second term describes RR. The run- 
away solutions LAD predicts can be avoided through 
asymptotic boundary conditions [T3], performing a re- 
duction of order to obtain the Landau Lifshitz (LL) equa- 
tion using an alternate equation or (since run- 
aways are nonperturbative in the classical electron radius 
To = e^/47rTO) simply by expanding the orbit in tq. 

In order to compare directly with QED, it is helpful to 
go back to the classical equations of motion and simply 
solve them perturbatively (rather than eliminating the 
radiation field to obtain LAD), in such a way that RR 
effects appear as corrections to the Lorentz force. To do 
so, define / := eF'^^^, let F be the dynamical electro- 
magnetic field with potential A, and j the current. The 
classical equations of motion are then, schematically |13) , 

i = /i + eFx , DA = ej . (2) 

(The external field obeys Maxwell's equations in vacuum, 
n^ext ~ 0). We solve ^ perturbatively in e, expanding 

X = xo + exi + . . . =^ j = jo + eji + ■ ■ ■ , , 

A = Aa + eAi + ... =^ F = Fq + eF^ + . . . ^ ' 

To zeroth order, we have DAq = and xq = f{xQ)xo, 
the Lorentz force equation. So, the particle moves ac- 
cording to the Lorentz force but, with appropriate initial 
conditions, Aq ^ Fq = and there is no radiation. To 
first order in e, we find a homogeneous equation for Xi. 
Since initial conditions can be fulfilled by xq, xi can be 
set to zero, and the particle's motion is unaffected. At 
this order we also have DAi = jo, and hence a nonzero 
radiation field -Fi sourced by xq, i.e. by a particle moving 
under the Lorentz force. Radiation is therefore created 
at order e, but field observables are typically of order 
e^, so the lowest order radiated energy, for example, is 
+ ^ e^jQ, which docs not contain recoil. 
At second order in e, one finds DA2 = F2 = and 

X2 = f{xQ)x2 + X2df{xo)xo + Fi{xa)xo . (4) 

This inhomogeous equation yields a nonzero X2- The par- 
ticle's orbit is corrected due to the term Fiio, i.e. by the 
fact that the particle has emitted the radiation Fi; this 
is radiation reaction. It appears at order in the elec- 
tron's motion, as expected. At third order, one finds 
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that F3 0, sourced by the radiating particle, i.e. X2 
contributes to the radiation field. Hence recoil effects 
appear in the radiated energy first at order e**, through 
the F1F3 ~ joj2 cross term. 

Example:- We illustrate the above for a background 
plane wave depending on invariant phase :— k.x with 
=0 [T]. We have fc^ — lju^ for a; an inverse-length 
scale, say a central frequency, and we can choose coordi- 
nates such that k.x = wx^, lightfront time. The trans- 
verse vector a'^ {k.a' — 0) is the normalised electric field, 
e£'j_ = uia'^. We consider pulses, so that E^{(j)) is ei- 
ther nonzero only in a finite 0-range or vanishes asymp- 
totically, but is otherwise arbitrary. The momentum 
TT^ := mig of a particle with initial momentum is, 
according to the Lorentz force [TH [TS] , 

f ,^ /,x 2p.a(0) — a^((/)) , , , 

and one has that k.p is conserved, hence cx t. (From 
here on, an integral without variable is over lightfront 
time d(j), tt = 7r((/)) and tt = 7r(oo) is the final momentum 
as given by the Lorentz force in j5|.) We proceed to • 
LAD, LL and the equations in pMoJ, while giving different 
results for X2 within the background, all agree on the 
first RR correction to the final momentum of a particle 
passing through the whole pulse; this is Stt^, with 

An explanation of why all the equations agree on this 
asymptotic result, to this order, is as follows. The first 
term in ^ can also be obtained by inserting the Lorentz 
orbit Xq into Larmor's formula for the total radiated mo- 
mentum. The second term in (|6|, representing momen- 
tum taken from the background, then follows by momen- 
tum conservation. 

Classical to quantum:- A direct, perturbative solu- 
tion of the classical theory has yielded the following. To 
zeroth order in e, a particle is accelerated by an external 
field but does not emit. To first order, the emitted radi- 
ation of the particle is accounted for. To second order, 
the effect of this emission on the electron is accounted for. 
The observables at this order include the corrected elec- 
tron orbit, and the emission spectrum of a particle mov- 
ing according to the Lorentz equation. At third order, the 
impact of recoil on the radiation field is accounted for, 
and enters the field's observables at order e*. 

Where will we find such effects in QED? Consider QED 
with a background field. In the Furry picture, interac- 
tions between quantised fields are treated perturbatively 
as normal, while the background is treated as a part of 
the 'free' theory and, classically, effects the Lorentz force 
on the fermions [T5HT7] . (The Furry picture therefore has 
the usual Feynman diagram expansion but with the free 
fermion propagator replaced by that in the background.) 



Hence, perturbative QED in the Furry picture is an ex- 
pansion in powers of a = e^/(47rfi,) beyond a 'free' the- 
ory without recoil, but with the Lorentz force. It is the 
quantum equivalent of the expansion used above. Thus, 
lowest order RR corrections should come from (the clas- 
sical limit of) order Furry picture processes in QED. 
We confirm this below. 

Radiation reaction from QED:- Consider the follow- 
ing simple scattering experiment. An electron is collided 
with a laser pulse, and its momentum is measured af- 
ter exiting the beam. The classical theory, using ([!]) or 
otherwise, predicts that the electron will have a certain 
momentum; -k + Stt from (|6| for a plane wave background. 
The quantum theory predicts that, repeating the exper- 
iment, a distribution of final momenta will be measured. 
The expectation value of this distribution can be calcu- 
lated in QED; it is the expectation value of the electron 
momentum operator in the final state describing the scat- 
tering process. Expectation values, unlike probabilities, 
have natural classical limits; this is why approaches based 
on them, such as the in-in formalism |18j , are common in 
studying how classical RR arises from QED [TH (THl HO] ■ 
Such approaches, though, are not as widely used in par- 
ticle physics as scattering amplitudes, so it is important 
to understand how RR arises in the latter. 

Returning to our experiment, we note that measure- 
ments are unlikely to be made within the fields of the 
laser, and so asymptotic results are sufficient. In this 
case, we can easily relate the expectation value of in- 
terest to familiar 5'-matrix elements and Feynman dia- 
grams, as follows. We begin with an incoming electron 
state evolve it in time through a background field, 
and calculate the expectation value of the electron mo- 
mentum operator, H^, in the evolved state. (In [21] . 
the position shift of a scattered particle was compared 
between quantum and classical theories. Here we com- 
pare momenta, which arc the natural variables in scatter- 
ing calculations.) If S is the S'-matrix, the expectation 
value is (H^) = {i\S^Ilf^S\i) . Inserting a complete set 
of asymptotic out states | / ) with H| / ) ^ pf\ f), we can 
write the expectation value in terms of S'-matrix elements 
Sf, = {f\S\i) as 

/ 

neatly relating the in-in and scattering formalisms. The 
sum in ([t]) runs over all Fock-sectors of free particle 
states, and can be interpreted as a sum of amplitudes 
for the processes e" — > e"+ anything. Each differential 
probability |S'/ip in ([t]) has a Furry- Feynman diagram 
expansion in powers of e. The zeroth order contribution 
to ([7]) comes from the one-electron sector, i.e. the process 
e^ — > e^ without emission, at tree level. Classically, it 
describes only the effect of the Lorentz force. The next 
contributions to (H^J come at order from two pro- 
cesses, see Fig. [l] First, photon emission at tree level. 
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FIG. 1. Feynman diagrams contributing to radiation reaction at lowest order. (A double line is the fermion propagator in a 
background field.) 



mod-squared. Second, the cross term of e~ — > e~ to one 
loop, describing the self-field of the electron. Together, 
these processes give a nonzero contribution to the final 
electron momentum and (H^); this is quantum radiation 
reaction to lowest order. The loop cannot be neglected. 
First, simply for consistency; one cannot keep one order 

contribution and drop another. Second, without the 
loop, (n^) would be infra-red divergent in general. It is 
a standard result that only the inclusive sum of the loop 
and emission processes is IR finite and contributes to ob- 
servables [22]. Such sums are automatically included in 
([t]); this is an advantage of considering (11) rather than 
individual diagrams. 

Example:- We now illustrate the above for a plane 
wave background, take H ^ and recover the classical re- 
sult As above, the zeroth order contribution to (11^) 
comes from tree-level scattering without emission, and 



yields (11^) 



[23], the Lorentz force result. We pro- 



ceed to order e^. Photon emission e"(p) — >• e^{p')+j{k') 
in a plane wave, called nonlinear Compton scattering, 
yields an outgoing momentum = p' with 



fc'.TT 



Pu 



k.{p~k') 



(8) 



and this should be inserted into the integrated probabil- 
ity of emission to obtain its contribution to (H^), see [l5] . 
Using similar methods, one finds that the loop in Fig. [T] 
leads to the same contribution as emission, except that 
the final momentum p^ — — tt rather than ([s]). The total 
contribution to the final electron momentum at order e^, 
call it (5(11^), is then (tTj = 7r((/)j), aj = a{(j)j)) 

"^dki f d^/cl k.p' 



4nh J k'_ J (27r)2 k.p 





i^fj. - Pp) \ d(?!>id(/): 



exp 



i f k' .TT 

Ti J "fcy 



and spin effects appear in 



9{t) 



1 



4(1 + t) 



where t 



k.k' 
k.p' 



(10) 



The result ^ is finite, nonzero and has support on the 
difference between the electron's momentum following 
photon emission, ([8|, and the Lorentz force (no recoil) 
result, TT |24j . It is due to combined photon emission and 



self-energy effects. This is quantum radiation reaction in 
a plane wave background. 

It is now known that momentum integrals such as those 
in ([9]) can be performed analytically and exactly [25| . 
opening the way for efficient calculations of cross sections 
in plane waves. Our objective here, though, is to take the 
classical limit of d{n^). To do so, we note that in a plane 
wave, it is the emitted photon's longitudinal momentum 
k.k' which is important. This 'breaks the symmetry' of 
motion associated with the Lorentz force [26] , namely the 
conservation of k.p. This is a sign of RR, as is seen explic- 
itly when one solves LL in a plane wave |27j . Hence, we 
will rewrite (|9| to highlight its dependence on k.k' . Let 
k'^ = (k.k' /k.p)r ^ and change variables k'^ — > r^, and 
k'_ ^ t — k.k' /k.p' , the ratio of longitudinal momenta in 
the final state. Finally, noting that photon momentum 
has no classical analogue, we rescale t ^ ht, which cor- 
responds to using wavenumber rather than momentum. 
This makes all dependencies on h manifest. The result is 
that Q is equal to (r_ = p_ since k' is on-shell) 



e 

4tt 



dt 



(1 + ht)^ J {2tt) 



d^r 



1 + ht 



r.TT 
k.p 



(11) 









l'd(t)id(j)2 exp 


t f ^ 


J k.p 





' — g{ht)[a2 ~ ai] 

r.TT2r.TTl 



where t is given by the following simple combination of 
final (scattered) momenta. 



ht = 



k.k' 



k.k' 



k.p' k.p — k.k' 



(12) 



(5(n^) depends on h only through the combination l-\-ht, 
so the classical limit h can be taken simply by ex- 
panding these factors in (111. From ([T2| , this is equiva- 



lent to assuming that k.k' <^ k.p' (implying k.k' <C k.p), 
i.e. that the momentum carried away by the photon is 
small compared to that of the electron. The classical 
limit therefore corresponds to neglecting quantum effects 
associated to emission of high energy photons. It remains 
to evaluate the t and integrals at h — 0. Noting that 
(5(n^) is real and symmetric in 0i o 02, the t integral 
can be symmetrised and immediately gives a delta func- 
tion for the two (/)-integrals, 



dicos t 



r.TT 
k.p 



k.p 



S{(t>2 - 01 



(13) 



4 



This result means that while quantum RR ^ is coher- 
ent in (lightfront) time, coming from an S'-matrix element 



squared, classical RR is incoherent; (13 1 is the statement 



that interference terms drop out and only the incoher- 
ent piece of the (/)-integrals in ^ survives as S — 0. 
We therefore see a simple aspect of decoherence, which 
is intrinsic to quantum-to-classical transitions [28j . Now, 
using ( 13 ) to eliminate one of the 0-integrals, the remain- 
ing Tj^ integrals become elementary, and we obtain 



a'H TT,, 



k.p 



Stt^ , (14) 



recovering the classical result ^ . We have now obtained 
classical radiation reaction directly from QED |29) . 

Discussion:- Related calculations appear in |30j . 
which partially recovers classical RR in a high energy 
approximation, and in |11| . which concludes that QED 
supports the classical equation in [5]. (We will address 
this elsewhere.) In perturbation theory, k.6n from ([6| 
can be obtained from ordinary Compton scattering if the 
incoming photon flux is equal to the plane wave's energy 
density a'^ [14] . That paper finds that quantum and 
LAD results differ for large recoil. This is explicit in our 
approach; accounting for larger recoil requires retaining 



higher order terms in the expansion of (111, and these 
terms are of higher order in H. 

We now turn to higher order RR effects. The sum 
in ([t]) is incoherent in particle number (as in |10j). but 
each process is coherent in time. Quantum RR thus re- 
tains some coherence, in general, rt-photon emission con- 
tributes to RR at order e^", but so do all other processes 
of the same order, e.g. pair production, loop corrections, 
counterterms and so on. These cannot be neglected in 
general. First, for consistency. Second, they remove UV 
and IR divergences. Third, unitarity can be violated 
without them. It was noted in [10] that photon emis- 
sion probabilities in plane waves can exceed unity. Such 
problems do not appear here; the terms which violate 
unitarity are removed by consistent inclusion of, e.g. the 
above loop, see also pS] . 

Further, our calculation shows that quantum processes 
can also play an important role in the classical limit. 
We find that the loop does not contribute directly to 
classical RR (which is natural) because it cancels against 
a term coming from the emission diagram; however, the 
cancelling terms are 0{l/h) and would otherwise diverge 
in the classical limit. Hence it seems that the loop is 
needed for both quantum and classical RR. See [3T] more 
on h and the classical limit. 

Finally, it is now clear (even from our classical discus- 
sion) why one-photon emission spectra do not exhibit RR 
[SI [9] ; at order , such effects are only visible in the elec- 
tron sector. The reason is that the probability of emitting 
a photon at all is already order , so one might say that 
a single emitted photon does not know that the electron 
it leaves behind has recoiled. RR effects appear in the 



p/ioton spectrum at higher order, following, e.g., multiple 
emissions. For the electron, there is no essential differ- 
ence between single and multi-photon emissions. Both 
contribute to RR. 

Conclusions:- Classical radiation reaction emerges 
directly from the h ^ limit of QED. We have identi- 
fied, in general, the contributing processes and diagrams; 
to lowest order, recoil effects in the electron momentum 
come from one-photon emission and one-loop self-energy. 

In an explicit example, we have seen how the recoil of 
an accelerated, emitting electron recovers classical radi- 
ation reaction in the low-energy limit. The advantage of 
using a plane wave as the background in this example is 
that no further approximations need be made; one sim- 
ply takes /i — > 0, thus eliminating potential ambiguities. 
Our calculation allows for a clear identification of recoil 
effects and unites the usual S'-matrix approach to QED 
with 'in-in' methods and the calculation of expectation 
values. 

Asymptotic results following from the S'-matrix, as 
used here, are sufficient for comparison with experiments 
aiming to observe radiation reaction, since scattering 
products will be measured far from interaction volumes. 
(Previous position shift calculations which find agree- 
ment between quantum and classical theories are also 
for 'asymptotic' times [H].) Nevertheless, the extension 
of our calculation to finite time and higher orders (it will 
be interesting to explore renormalisation effects here) can 
distinguish between, and therefore rule out, different clas- 
sical equations which predict different motion within the 
background field. Our results will appear soon. 
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